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INTRODUCTION 


This  report  is  concerned  with  the  absorption  of  radiation  at  optical  wavelengths  i rji 
a  plasma  aerosol  ( [  1  ]- [3]).  We  shall  take  this  to  include  wavelengths  from  3  cm  to  0.3  fi 
The  term  aerosol  refers  to  a  gas-particle  mixture  such  as  that  found  in  the  atmosphere. 
Furthermore,  the  plasma  is  assumed  to  be  of  the  type  achieved  by  laser-induced  air 
breakdown  described  in  jl]  and  |3j.  A  number  of  conclusions  may  be  reached  just  by  an 
examination  of  the  absorption  mean  free  path  lu  (the  path  length  necessary  for  th+e^  inten¬ 
sity  of  the  impinging  radiation,  e.g.,  a  laser  beam,  to  fall  to  1  e  of  its  value).  More 
detailed  information,  such  as  the  potential  consequences  of  the  dissipation  of  absorbed 
energy,  depends  on  specifying  explicitly  the  geometry  (for  both  the  beam  and  the 
plasma),  characteristic  times,  and  the  mechanisms  of  plasma  containment  'maintenance. 


ABSORPTION  BEHAVIOR 

For  both  historical  and  practical  reasons  there  exists  a  plethora  of  formulae  for  the 
absorption  coefficient  /(,,  The  various  derivations  may  be  roughly  classified  into  three 
categories  those  using  the  Vlasov  transport  equation  (classical),  those  using  classical 
biemsst  rahliing  coupled  with  the  principle  of  detailed  balancing  (sennclassical).  and 
quantum  mechanical  In  all  except  the  least  rigorous  derivations,  these  techniques  result 
in  expressions  involving  unpleasant  integrals  Consequently,  approximations  are  made 
which  result  in  a  large  set  of  formulae,  each  valid  under  particular  conditions  appropri¬ 
ate  to  its  derivation  The  appendix  of  this  report  contains  an  extensive  summary  of 
those  which  are  pertinent  to  the  present  study 

It  can  be  seen  from  those  equations  that  the  temperature  dependence  of  /i^  is  quite 
('implex  Nevertheless,  in  considering  an  aerosol  plasma,  we  shall  find  ourselves  able  to 
pick  a  nominal  temperature  on  which  to  base  our  conclusions  We  first  note  that 
because  of  the  exponential  dependence  of  the  electron  partial  pressure  on  temperature  in 
equations  IA-22)  and  (A-30).  there  is  a  steep  rise  in  absorption  just  prior  to  complete  sin¬ 
gle  loni/aj^n  For  atmospheric  pressure  tills  occurs  at  about  T  1  5  eV  or 
17  ">()()  "K  At  that  point  the  electron  pressure  is  about  one-half  the  total  pressure, 
and  /i  _  is  given  by  (A-l!))  with  Z  I  It  can  be  seen  from  that  equation  that  /i  v 
decreases  with  a  further  increase  of  temperature  Tins  holds  true  until  about  'Jo. 000  "lx. 
when  double  loiii/atioii  begins  to  dominate,  and  //  .  begins  to  rise  again  (  1  .  p  213) 

I  a r i her  computations  require  the  entire  set  of  equations  (A-21 ),  (A-22).  (A-27),  and  (A- 
2M  until  the  double  ionization  is  close  to  complete,  at  which  point  ( A- 1 0 )  is  again  a  good 
approximation  but  with  Z  2  As  the  temperature  rises,  this  procedure  repeats  itself 
and  /<  passes  through  a  series  of  maxima  and  minima  The  ratios  of  adjacent  maxima  to 

minima  are  only  on  the  order  of  two  (cf  .  I  .  p  213)  so  that  it  is  not  unreasonable  to 

examine  only  the  maximum  values  Also,  as  noted  above,  temperatures  below  ionization 
m  iv  be  dismissed  as  having  insignificant  absorption  Finally,  at  or  near  breakdown  for 
aerosol  plasmas  the  temperature  is  in  the  range  of  1  to  2  eY  (}3j,  p  966),  so  that  as  a 

t  Tli.It  i-  10  <;||z  to  10"  <;Hz  "I  I  0  X  10  *  eV  I.-.  I  O  eV 
ft  /  i'  I  /i  whi'ti'  il  is  tin'  ali-'i pt I- -it  ( ■  .efficient 

ft*  s. .•  jn.it i  n  i  A  Jhl  I  -  i  l.iign  |'i i"~s hi i".  it  i in  i I'.ises  ("i  example  at  p  ,  ---  100,  total  10m- 

.■  iti  a  o  in  ,e  .<!  ■  in  T  2K. 000  "lx 


general  guide  to  absorption  we  may  assume  that  there  is  complete  single  ionization;  i.e 
we  take  T  =  17,500  °K  and  Z  —  1. 

Let  us  now  compute  u>p  and  i/m  in  order  to  determine  which  approximations  are 
valid.  The  ion  density  is  given  by 


n 


P 

2kT 


1.013  x  10s  patm 
2  (1.38  x  10“16)  (1.75  x  104) 


=  2.10  x  1017  patm. 


(1) 


Then,  at  1  atm 


u.’p  =  (4tre2  2.10  x  1017/m)1,2  =  2.58  x  1013  sec'1. 


This  is  equivalent  to  a  wavelength  of 


6.28  (3.0  x  IQ10) 
2.58  x  1013 


7.3  x  10  3  cm 


C2) 


(3) 


Also,  from  (A-9)  we  have 


v 


m 


diL  f  1 1  e*  (2  10  x  1Ql7) 

3  3  J  (k  1  75  x  104)3^2\/m 


~  (6.0  x  10n)  g  sec  ’. 


(■*) 


Thus,  <v  corresponds  to  X  «  -  0.3  cm  or  freq  100  GHz. 

Note  that  for  wavelengths  greater  than  0.01  cm  (i.e  .  frequencies  below  3000  GHz) 
there  is  total  reflection  of  impinging  radiation.  Furthermore,  since  «  vtp,  we  see 
that  for  frequencies  much  greater  than  3000  GHz  conditions  (ii)  and  (iii)  of  t  he  Ap  pendix 
are  both  satisfied.  Finally,  down  to  wavelengths  of  about  3.0  x  10  4  cm.  condition  (i')  is 
also  valid.  (More  precisely.  W  <>,  kT  for  X  ■  •  8.2  x  10  ’em  0.82  //,  and  tiw  •  I] 
for  X  »  8.5  x  10  6  cm  —  0.085  ft  )  In  the  region  from  3  x  10  s  cm  to  3  x  10  4  cm  we  use 
the  Kramers-Unsold  formula  (A-30b),  which  is  justified  since  h lj.  At  longer 
wavelengths  (i)  is  valid,  and  we  can  use  the  simpler  equation  (A-19e),  which  becomes 

/i„.  :  3.8  x  105g  p2tm  X-’  (5) 

To  evaluate  the  Gaunt  factor  g,  we  first  determine  Ac.  From  (A-3) 

A(-  3.8  x  103  X  (6) 


•) 


The  low-frequency  region  occurs  for  X  »  2.6  x  lCT^cm,  which  we  arbitrarily  set  at 
X>3xl0_3cm;  the  remainder  of  our  spectrum,  3  x  l(T3cm  to  3  x  10'7  cm,  falls  in  a 
region  for  which  we  have  no  simple  formula  for  g  (or  £0(oj)  ). 

Figure  1  contains  a  plot  of  lu  versus  wavelength.  In  the  figure,  we  have  marked  the 
region  where  the  Kramers-Unsold  formula  was  used,  and  the  factor  £o(w)  was  simply  set 
equal  to  two.  This  matches  the  value  obtained  by  (5)  at  X  =  1CT3  cm,  a  choice  of  £0  that 
differs  slightly  from  [1],  page  242. T  Observe  the  leveling  of  the  curve  at  X  equal  to 
2  x  10  s  cm.  At  even  smaller  wavelengths  (which  happen  to  be  outside  our  region  of  con¬ 
cern)  it  has  a  sawtooth  shape,  alternately  increasing  and  decreasing  ([4],  p  274).  The 
high-  and  low-frequency  regions  used  to  compute  g  in  equation  (5)  are  also  noted  in  the 
figure.  We  have  used  the  table  in  [4],  page  255,  to  compute  g  in  the  high-frequency 
region.  It  is  only  tabulated  down  to  ticj/kT  =  0.1;  i.e.,  up  to  X  =  8  x  10  4  cm.  From 
this  wavelength  up  to  the  low-frequency  region,  g  was  set  equal  to  2.0,  the  bottom  entry 
in  the  table.  This  procedure  seems  to  have  produced  a  good  interpolation  between  the 
tw'O  regions,  at  least  on  the  log-log  plot.  Even  though  the  table  and  (A-3)  include  the 
frequency  dependence  of  g,  the  slope  of  the  curve  in  the  figure  is  essentially  a  constant 
within  the  range  of  equation  (5).  The  accuracy  is,  of  course,  not  very  good  close  to  the 
plasma  wavelength  Xp  —  7  x  l(T3cm. 

Finally,  if  we  approximate  (5)  by  setting  g  =  2.0  (a  reasonable  value  for  the  optical 
region),  we  obtain  the  particularly  simple  formula 

/  =  8  x  1(T5 

2  \  2 
Palm  " 

which  is  essentially  the  same  as  that  in  3  Note  that  if  the  impinging  beam  is  not 
parallel  (e  g.,  a  focused  laser),  the  beam  divergence  should  be  taken  into  account  in  com¬ 
puting  absorption  distance. 

CONCLUSION 

We  have  seen  that  at  atmospheric  pressure,  there  is  significant  absorption  only  at 
wavelengths  above  3  x  lO^1  cm  —  3  //  (Also,  the  medium  is  essentially  a  total  reflector 
for  X  >  10  “cm  100  /<  )  One  obvious  way  to  extend  this  range  is  to  increase  the 
plasma  pressure.  According  to  3;,  for  short  periods  of  time  during  breakdown  the  pres¬ 
sure  at  a  breakdown  site  may  lie  as  high  as  100  atm.  Whether  it  is  feasible  to  maintain 
a  pressure  of  this  magnitude  over  longer  periods  and  throughout  the  plasma  depends  on 
the  details  of  the  plasma  maintenance  and  containment  (Most  treatments  of  the 
maintenance  of  a  laser-induced  air  plasma  assume  close  to  atmospheric  pressure;  c.f.,  2], 
'5.) 

|  Raizer  sets  g  -=  2  5  for  a  C02  laser  (X  lO^cm).  and  he  sets  {-=07  for  a  neodymium  laser 
(X  ss  10  *  cm)  He  also  includes  an  "edge"  factor  of  1  2  to  1  5 

tf  For  X  smaller  than  the  low-frequency  region,  that  is,  for  X  from  3x  lO^cm  to  3  x  10"®  cm, 
ti^'.'kT  ranges  from  2  7  to  0  027 


bO  4C 


NOTATION 


c  =  2.998  1010  cm/sec  (speed  of  light) 

e  —  4.803  10"10  esu  (charge  of  electron) 

g,  gA,  gB,  etc.  =  Gaunt  factor  (cf.  equation  (A-3)  ) 

”1  =  statistical  weight  of  ground  state  of  ion  of  charge  m 
=  6.626  10_27/2 tr  erg-sec  (Planck’s  constant/2;r) 

=  1.055  1027  erg-sec 
I  =  ionization  potential  of  plasma  atoms 
=  Z'IB,  where  IH  is  that  of  a  hydrogen  atom 
k  =  1.380  10  16  joule/ °K  (Boltzmann’s  constant) 

{j  —  absorption  mean  free  path  (  =  1  //iu,) 
m  =  9.109  1CT28  grams  (electron  rest  mass) 
n.  =  number  density  of  electrons 
nm  =  number  density  of  ions  of  charge  m 
Note  that  Z  =  m  +  1 
n  =  total  number  density  of  ions 
n_  =  number  density  of  ions  when  ionization  is  total 
(i.e.,  assuming  all  ions  have  the  same  charge) 

P»(m  —  pressure  in  atmospheres 
p  —  pressure  [dynes/cm2] 

=  1013  106  x  Palm 
T  =  temperature,  °K 

u,m)  —  partition  function  for  mth  ion,  normalized  by  the  ground  state  energy 
x.  =  p./p.  where  p  is  electron  partial  pressure 
Z  charge  of  nucleus  for  hydrogen-like  atom 
"7  =  1.781  (ln(7)  =  Euler’s  constant) 

X  =  wavelength  in  cm 

=  total  absorption  coefficient  for  monochromatic  radiation 
of  frequency 

um  =  effective  collision  frequency  (for  momentum  c  hange) 
i>m  =  velocity-independent  approximation  of  vm 
u>  —  frequency,  radians/sec  =  2tr  frequency 
Note  that  u >  (sec'1)  =  1.518  1015  x  (tiw)ev 
=  plasma  frequency  in  radians/sec 
=  (4n-e2n./m)l/2 


G 


APPENDIX:  ABSORPTION  COEFFICIENT  FORMULAE 


FREE-FREE  TRANSITIONS 

One  of  the  more  useful  expressions,  from  a  viewpoint  of  generality  and  eompu t  abil¬ 
ity,  for  the  absorption  coefficient  of  a  plasma  is  ([4],  p  259-260) 


ku  ( 1  e 


w/kT)  [cm  1 


1  f)7r_ 


7  2 

1  /  c  n  .  n 


( mkT)  •  me  ti  . 


I  A  1  l 


I  A- 2) 


and  g  is  a  so-called  Oaunt  factor  given  by  (  I  p  255.  6  p  488-  191,  7  equation  (5-591) 


x/a 

- In  A(. 

for  T 

(8.9  x  10’)  /.-  and  Ar  ■  1 

(A  dal 

8  -  • 

0(1) 

for  T 

<  (8.9  x  1 0° )  7?  and  Ac  =  1 

( A-3b) 

n/3  . 

— —Ill  Ay 

for  T 

.">  (8.9  x  10*)  7r 

(A  3c) 

where  |cf  ecpiations  |  A  - 1 1  )  and  (A- 12)  ) 


A( 

A,, 


2  1 2kT) 


a 


•a  'J 


•Ml  x  10 


2  9  1  x  1 0‘ 


T 


*;■’  ~Ze* v  III.. 

IkT 

1 1 

in  a  a  1 1'  'ii  ( )(  1  I  is  in  be  in  t  ejprei  ed  as  "  <  4  the  order  I 


h.ipiat  lolls  I  A-  I  I  and  (  \-2)  ale  derived  ill  I  siiii  K  lass|<  allv  ,  I  e  ,  brellisst  la  ll  III  II  g 
radiation  for  an  electron  in  a  coulomb  potential  is  computed  classically,  the  coefficient 
x  deduced  from  an  application  of  the  principle  of  detailed  balancing  plus  Planck's  for¬ 
mula  for  radiative  equilibrium,  and  the  factor  in  (A-l)  added  to  correct  for  stimulated 
emission  |  he  assumptions  used  in  the  derivation  are 

( I )  I  I  ee-l  I  ee  t  rail  si  I  |o||s  are  t  he  dolll  I II  a  ll  t  absorpt  ion  Itiech  all  Is  III 
(ii)  -c  •  •  -Cj,  where  is  the  plasma  frequency. 

(in)  ss.'  •  ■  u,n  where  //m  is  the  transport  collision  frequency. 

♦The  ( ;  ,i  1 1  ii  I  f.n  I.  i  w  hit  li  i-  n  I  .jei  n  eii  mill  is  <|e,  iis.-etl  hr|.  .w 


(iv)  hot  «  kT. 

Note,  there  is  not  a  simple  expression  for  g  which  is  accurate  in  the  crossover  region 
Ac  1 . 

Condition  (i)  is  tantamount  to  ignoring  energy  absorbed  in  photoionization  of  the 
atoms  or  ions  (bound-free  transitions).  It  is  intuitively  clear  that  a  necessary  condition 
for  this  approximation  is  that  the  photon  energy  1w  be  small  in  comparison  with  I,  the 
ionization  potential  of  the  atoms.  Furthermore,  it  can  be  shown  ([4]  p  271-272)  that, 
since  under  these  conditions  photons  are  only  absorbed  by  highly  excited  atoms,  the 
ratio  of  the  contribution  from  bound-free  transitions  to  free-free  is  e1'u'^kT-l  ([4]  p  27 1  - 
272).  Thus,  a  sufficient  condition  for  (i)  is 


(i'a)  hu.’  «  I 

(i'b)  tico  «kT 


Linder  these  circumstances  (A-l)  reduces  to 


Pu  =■ 


1 67T" 

3 


2?r  *  "  e6Z’*n+n 

c(mkT)3/2J* 


(A- 4) 
(A-5) 


(A-6) 


which  is  equation  (6.1)  of  [l]  (often  referred  to  as  Kramers  formula).  On  the  other  hand, 
if  ttw>  kT,  the  semiclassical  treatment  of  free-free  transitions  and  even  the  (quantum) 
approximation  (A-3c)  lose  their  validity;  moreover,  in  that  case  bound-free  absorption 
dominates  ( [4]  p  272). 

Condition  (ii)  plays  two  roles  in  the  derivation  of  (A-2).  The  first,  which  is  rela¬ 
tively  transparent,  is  that  it  insures  that  the  index  of  refraction  is  close  to  one  (cf.  1  p. 
28,  [6)).  The  application  of  this  assumption  in  the  derivation  is  fairly  clear  It  enters 
into  the  energy  balancing  computation  via  the  relationship  between  the  energy  flux  and 
energy  density  ([1]  p  28),  which  assumes  that  the  speed  of  energy  propagation  is  equal  to 
the  speed  of  light  in  a  vacuum.  (Compare  [6],  equations  (12.17)  -  (12.18),  where  c  is 
replaced  by  the  group  velocity  and  n  is  not  one.)  Classically,  for  *.<  <  u;p  there  is  total 
reflection.  Formulae  for  the  absorption  coefficient  when  vC  is  in  the  neighborhood  of 
may  be  found  in  [8],  [9],  [1]  p  308,  and  [6 j  p  66-68.  The  less  obvious  role  for  condition 
(ii)  relates  to  Debye  shielding.  Formulae  (5.8)  and  (5.9)  of  [4]  are  obtained  by  integrating 
over  the  impact  parameter  b  for  electron  scattering  with  an  upper  limit  of  infinity.  This 
is  equivalent  to  ignoring  the  shielding,  and  it  is  valid  provided  u;  »  u;p  (see  ( 10] ).  We 
note  that  [8]  includes  shielding  effects  in  the  transport  equation,  but  only  gets  simple 
expressions  for  the  cases  uj  »  utp  and  w’  «  u>p. 


t  Heunstically,  we  note  that  contributions  at  large  b  only  affect  frequencies  u>  <  v0/b  (where  v0  = 
initial  electron  velocity)  If  we  take  the  cutoff  to  be  the  Debye  length  XD,  then  the  error  produced 
by  extending  our  integration  to  values  of  b  larger  than  XD  is  negligible  at  those  frequencies  satisfy¬ 
ing  w  »  v0/XD  (kT/m)1/2/XD  —  wp 
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Inequality  (iii)  expresses  conditions  under  which  the  correlation  of  successive  colli¬ 
sions  of  an  electron  with  ions  and  its  attendant  interference  of  bremsstrahlung  radiation 
may  be  ignored  ([*4]  p  258]  The  quantity  i>m  is  the  collision  frequency  for  the  momen¬ 
tum  change;  i.e.,  the  mean  change  in  momentum  per  unit  time  due  to  collisions  is  miq„v 
([l],[6j).  In  general  vm  is  a  function  of  speed  v  (cf.,  [6]  p  51-52);  however,  for  the  pur¬ 
poses  of  (iii)  we  may  use  the  expression  (A-9).  When  uj  is  of  the  order  of  vm,  an 
appropriate  formula  is  (A-7),  which  includes  the  effects  of  correlated  collisions;  this 
reduces  to  the  simpler  expression  (A-8)  provided  um  is  essentially  a  constant  over  the 
range  of  velocities  which  are  significant  in  the  integral. 

(iv)  The  derivation  of  (A-3a)  and  (A-3b)  involves  two  integrations.  The  expression 

for  bremsstrahlung  is  obtained  by  an  integration  over  the  electron’s  impact  parameter 

11;.  balancing  leads  to  an  absorption  coefficient  which  is  a  function  of  electron  initial 
velocity  v0,  and  this  expression  is  integrated  over  the  Maxwell  distribution  to  obtain  kuj. 
The  bremsstrahlung  formula  is  simplified  via  two  approximations,  one  at  low  frequencies 
(  •I  equation  (5.9))  and  one  at  high  frequencies  (  -4  equation  (5.8)).  Upon  balancing  and 
integration  over  v0  (cf.  (>!),  these  two  approximations  lead  to  (A-3a)  and  to  g  =  1 , 
respectively.  However,  the  conditions  of  approximation  are  velocity  dependent,  so  that 
one  must  show  that  the  contribution  to  the  integration  over  v0  outside  the  region  of 
approximation  is  insignificant  Oster  10  does  this  for  low  frequencies  and  shows  that 

the  condition  lor  validity  is  A(.  •  •  1  Ginsberg  (  f>  p  189)  seems  to  feel  that  the  pro¬ 

cedure  is  justified  in  both  cases.  (Actually,  4  only  treats  the  high-frequency  case,  which 
leads  to  the  Kramers  formula  (A-2)  with  g  1.)  We  note  that  even  if  we  accept  the 

value  g  =5=  I  in  the  high-frequency  region  where  Ac  •  <  1  (but  still  requiring  ti  a.'  «  1,  of 
course),  there  is  still  a  gap  for  frequencies  where  At  ^  1 .  A  rough  table  of  Gaunt  factors 
appears  in  -4  ,  p  255.  The  two  other  approximations,  involving  T,  separate  the  classical 
region  from  that  in  which  a  quantum-mechanical  derivation  is  necessary.  Quantum- 
mechanical  effects  become  important  where  the  (correspondence  principle)  quantum 
minimum  impact  parameter  ti  mv„  approaches  the  classical  minimum  impact  parameter 
Xe*  mv,f  (c.l.  9  )  In  the  quantum-mechanical  treatment,  the  impact  parameter  does 
not  appear  explicitly;  however,  the  results  for  a  single  electron  interaction  do  have  to  be 
integrated  over  the  velocity  vH  The  result  may  be  approximated  bv  the  simple  expres¬ 
sion  ( A-3c)  under  the  condition  Aq  •  \^Af;  i.e..  T  •  5x  10  s  X"  10  Oster  10  also 

shows  that  at  low  frequencies  (i.e.,  Ar  1)  the  general  quantum-mechanical  expression 
reduces  to  the  classical  formula  (A-3a)  under  the  reverse  inequality  T  5x  10’  X" 

Both  these  approximations  also  involve  the  condition  h  _  kT  The  region  in  which  a 
classical  approach  is  permissible  may  also  be  obtained  heurist teallv  by  subsjjtuting  an 
average  velocity  ■  v„  (kT  m)1  "  in  the  relation  Xe"  mv0*  h  mv0  Observe 

that  in  (A-3)  we  have  interpolated  over  the  crossover  region,  ic  .  set  the  switching  point 
at  Ac  -  Aq,  which  gives  T  8.9  x  10”  X"  Oster  seems  to  feel  that  this  interpolation  is 
justified 


t  Hv  wlnih  IS  ine.int  all  aver.ige  .  ve|  the  iiii|m<  I  parameter  wiOi  a  fixed  initial  ve|. ..  itv  •  f  the  rlo  ■ 
10  .|| 

tt  An  exariniial  ion  of  the  denvati.  n  "f  (A  -  d<  )  mdn  ates  that  a  natui.il  <h.>ne  for  v0  -  is 
I -’kT  im)1  2  [  10)  This  leads  to  the  inequality  T  •  -  8  x  10* 
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EXPRESSIONS  FOR  um  AND  g 

A  classical  treatment  of  the  problem  utilizing  the  transport  equations  yields  ([1)  p 

42,  [6]  p  51-52) 


Vv  - 


47re* 


3mco;‘  o  w+tC 


V  V3— (0v'2)dv 


t 


dv 


(A-7) 


which  is  valid  under  the  condition  (i').  The  factor  ur/(oj2  -f  v£)  reflects  the  effects  of 
correlated  collisions  ( [l]  p  29-34,  [4]  p  257);  i.e.,  condition  (iii)  is  no  longer  necessary.  If 
um  may  be  treated  as  a  constant  i>m,  then  (A-7)  can  be  integrated  by  parts  to  yield 


4^e‘n_  L>m 

m°  u,’"  -t-  i 


(AS) 


Alternatively,  we  may  define  i/m  as  that  value  for  which  (A-8)  holds  in  the  limit 
uJ  »  Pm'  Applying  this  definition  to  (A-7)  with  <t>  the  Maxwell  distribution  yields  (  6  p 
52-53' '  and  p  490) 


(\  1/2  rjZ  4 
111  z  e  11  ^ 

3  J  (kT)3  2\/m 


( A-9 ) 


Dropping  i>2  in  the  denominator  of  (A-8)  and  substituting  (A-9)  in  that  equation,  we 
arrive  at  (A-6). 

The  precise  expression  for  g  bears  some  discussion.  In  elementary  treatments  (cf. 
;lj  and  [4')  it  is  simply  unity  A  very  coarse  analysis,  employing  only  average  quantities, 
may  be  performed  as  follows.  Let  v=(kT/tn)1  *  be  the  average  electron  speed.  Let  r,.. 
the  coulomb  radius,  be  the  minimum  radius  to  which  an  electron  of  average  energy  can 
approach  and  still  escape  capture;  i.e.,  Ze*  r,.  —  mv’  2  kT  2.  Finally,  let  the  coulomb 
cross  section  bi  ctc  =  7rrc*.  and  i>m  —  n  +  v<r(,.  It  then  follows  that 


47m,  Z‘V 

"  ( kT)3  2v/m 

i.e., 

gA  -  3(3/2»r),,a  =-  2.1  .  (A- 10) 


A  more  precise  treatment,  employing  transport  equations  including  the  Boltzmann 
collision  term,  but  nevertheless  averaging  tz|n  over  the  velocity  distribution,  yields  the 
expression  (A-9)  for  i/m,  with  g  given  by  (jl2j  p  251  and  348,  [13j  p  258) 

t  #v)dv  is  the  electron  density  for  speed  v 

tt  Note  that,  in  equation  (6  13)  of  [6|,  v  —  (8kT/m),/2  (p  52)  and  ln(0  37kT  )  (2  i)  ln(  )s 

—  ( 2/3)( w/x/3)  g  which  results  in  (A-9) 
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gp  = 


3  (kT)3/2 

v/Tir  Ze3n.1//2 


(A- 11) 


y/3 


in 


Note  that  this  A  differs  from  that  of  [4|  on  page  -118  by  a  factor  of  one-ha!f. 

Alternatively,  deducing  um  from  a  comparison  of  (A-7)  and  (A-8)  in  the  limit 
u,'»i'rn  yields  a  Gaunt  factor  with  a  frequency  dependence  ([6]  p  490,  [9]).  In  particular, 
we  have  (cf.  [6)  p  490) 


Sc 


_  vi,n  2(2kT); 


3/2 


tr  75/2Ze-vV 


uJ 


(A- 1 2) 


where  7  —  e0j"  --  1.781.  The  same  expression  for  g  is  found  in  [7 j  (p  148)  and  jlO'j. 
We  note  that  these  formulae  are  essentially  classical  and  apply  only  to  temperatures  in 
the  range  ([6!  p  488,  [9j,  [10]) 


2  71'  j 

i  «.  - - - —  =5=  5  x  10°  for  Z  —  1. 

2  h  Jk 


(A- 13) 


In  the  opposite  extreme,  at  high  temperatures,  we  have  ([6j  p  491,  [9]^,  [10),  [1 4j  p  246) 

\/3  i  4kT  ,A  . 

gD  =  - In  — —  .  (A- 14) 


Vh, 


which  must  be  derived  from  quantum-mechanical  considerations.  All  of  these  expres¬ 
sions  (  (A- 11)  through  (A-14)  )  are  only  valid  for  frequencies  in  the  region  in  which  the 
argument  of  the  logarithm  is  much  greater  than  one. 

The  variations  in  the  equations  for  g  tend  to  be  the  result  of  the  somewhat  arbi¬ 
trary  cutoffs  involved  in  the  classical  approximations;  nevertheless,  their  role  can  provide 
physical  insight.  For  example,  [9]  provides  the  less  precise  expression 


gE 


v/3  ,  vtTT 

- In  - —  , 

tf  Ze"u > 


(A- 15) 


which  reduces  to  an  expression  very  close  to  (A-12)  if  the  thermal  velocity  vT  is  taken  as 
( kT  m)  *.  Another  approximation  is  to  avoid  the  frequency  dependence  by  replacing  u! 
with  -cp  in  (A-12).  This  leads  to 

t  This  ran  be  written  7  =  e1  ,  where  7  is  Euler’s  constant  Care  should  be  taken  inasmuch  as 
some  references  (cf  ,  [7]  )  use  the  notation  7  for  7* 

ft  |9j  finds  A<3=vi{mkT)l/2/tiw,  which  requires  a  value  for  the  electron  thermal  velocity  vT  This 
differs  from  (A-16)  by  4/1  781  =  2  25  in  the  argument  of  the  logarithm  if  vT—  (kT/m)1/!! 


II 


(A- 16) 


73  2(2kT)3/2  vW 

"  q5/2Ze27^  (4«2n.)I/2 


73 


In  .38 


(kT)3/2 
Ze3n  >/* 


If  ~j  is  set  equal  to  one,  then  this  matches  (A-ll)  up  to  replacement  of  the  factor  3/ TT/t 
0.85  by  0.80  .  We  have  adopted  (A-12)  in  conjunction  with  (A- 1 )  as  the  most 
rigorously  derived  under  conditions  (i)  -  (iv).  Note  that  (A-16)  does  not  reduce  to  the 
Gaunt  factor  given  by  Raizer  ([l]  equation  (6.1))  even  though  the  footnote  to  that  equa¬ 
tion  might  seem  to  indicate  that  should  be  the  case.  In  fact,  if  we  multiply  by  2/3  and 
introduce  the  exponent  3^2  inside  the  logarithm,  Raizer’s  expression  becomes 


Sc; 


73  2 

n  3 


In  8 


(kT)3/2 


(A- 17) 


which  differs  from  the  various  expressions  in  the  exponent  of  Z  as  well  as  the  numerical 
coefficients.  We  do  not  know  whether  this  expression  of  Raizer  is  erroneous;  however,  it 
does  yield  reasonable  results  in  the  region  of  application.  (If  T  —  15,000  °K  and  p  = 
atmospheric  pressure,  then  the  electron  partial  pressure  pe  is  about  106/2  dynes/cm2,  n_ 
—  Pe/kT  =  2.42  x  1 017,  and  we  have  g  =  2.23  for  (A-17)  while  g  =  1.67  from  equation 
(A-16).)  Most  likely  (A-17)  represents  an  approximation  to  the  Gaunt  factor  for  optical 
frequencies  inasmuch  as,  under  those  circumstances,  Ac  1,  and  (A-12)  is  unacceptable. 

IMPLICATIONS  OF  THE  SAHA  EQUATION 

To  actually  compute  /<w,  we  need,  in  addition  to  values  for  the  temperature  and  fre¬ 
quency,  the  ion  and  electron  densities  n+  and  n_.  Both  of  these  may  be  computed  from 
the  pressure  under  the  appropriate  thermodynamic  assumptions.  For  a  perfect  gas  in 
equilibrium  (i.e..  T+  =  TJ,  we  have 


P+  =  n  +  kT 

p_  =  n.kT  (A- 18) 

p  =  p+  +  p_  +  Pn  • 


where  p+,  p_,  and  pn  are  the  partial  pressures  of  the  ions,  electrons  and  neutrals.  We 
still  have  too  many  unknowns;  however,  if  we  assume  total  ionization,  then  pn  =  0, 
n_  =  Zn  +  ,  and  (A-6)  becomes 


16rr2 

2n 

1/2  K  v  2 

eBg  ZP 

~  3 

3 

m3/«c(kT)7« 

16tt2 

2tt 

'n  eV  ZJ  p7 

3 

3 

m3/2c(kT)7/2  (l  +  Z)2  w2 

(A- 19a) 


(A- 19b) 
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3  77  x  1029 


Z3 


Patri 


=  0.16 


(1  4  Z)2  (T/104)7 
Z3  _ g 


2 

Pain 


(1  4  Z)2  (T/104)7'8  (tw)e2 

-  z3 


1.07  x  10‘ 


(1  4  Z)*  (T/10 


(A- 19c) 
(A-19d) 
(A-lfle) 


More  generally,  a  gas  will  exist  in  various  stages  of  ionization.  The  particle  densi¬ 
ties  for  a  monatomic  gas  may  be  described  by  the  Saha  equation  (!4j  p  192-207) 


■'  in  1 1 n 


jlni-t  1  j 
.>1) 


mkT 
2ir  h  2 


3  -2 


Iim  l  tT 


(A- 20) 


where  nm  is  the  number  density  of  ions  of  charge  m,  m  is  the  number  density  of  elec¬ 
trons,  urn  is  the  partition  function  for  the  mth  ion  normalized  by  the  ground  state  energy 
(see  j-11  equation  (3.-43)),  and  lrn  is  the  ionization  energy  to  remove  the  mth  electron  after 
the  first  m-1  have  been  removed.  The  particle  densities  may  be  related  to  the  partial 
pressures  and  thence  to  the  total  pressure;  however,  solving  the  resulting  set  of  equations 
even  for  the  case  of  a  gas  consisting  of  a  single  element  can  be  quite  involved  More 
specifically,  in  addition  to  (A-20),  we  have 


nrmx 

11  Z>n.  4n  (A-21a) 

0 

nmiLX 

n.  -  £mnm-  ( A  -  2 1  b ) 

0 

where  n  is  the  total  number  density  of  particles,  (liven  n,  this  provides  nmax  -+-  2  non¬ 
linear  equations  in  the  same  number  of  unknowns. 

Things  simplify  considerably  when  only  one  level  of  ionization  is  present,  i.e  .  when 
nm  0  except  for  one  set  of  adjacent  values,  m  -  mn  and  m  mn  +  1  (Note  that 
Z  m„  4-  1.)  The  simplest,  and  for  our  purposes  most  useful,  case  is  when  Z  1  Then 
m0  —  0,  so  that  n  +  ,  which  by  definition  is  nmo+1,  equals  n  ,  and  no  =  n  2n  .  In  terms 
of  partial  pressures  the  left-hand  side  of  (A-20)  may  be  written 

IM>  P  ~ 

kT  p0  kT  (p-2p  ) 

furthermore,  if  we  let  x_  --  p./p  and  approximate  the  partition  functions  by  the  ground 
t  A  brief  treatment  of  multiple  ionization  nuv  tie  f..und  111  [4, 
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state  statistical  weights  g*m*.  we  get  equation  (6.-17)  of  1 


g(11f  m 

)3  *  (kT)5/2  -t./KT 

g*0)  ^  2?rh  2 

I  p 

=  6.7  x  103  C 
g(0! 

(T,  104)V2e  kT 

(A- 22 

Patm 

These  expressions  may  be  combined  with  (A-6)  to  give  a  formula  for  the  absorption 
coefficient  which  (assuming  conditions  (i').  (li).  (iii).  (iv)  )  is  also  valid  for  weakly  ionized 
gases: 


1.07  x  10' 


(T  101 


Patm  ^ 


(A-23) 


Because  of  the  exponential  in  (A-22).  at  some  temperature  depending  on  I } .  the  gas 
will  experience  a  sudden  rise  in  the  electron  partial  pressure  to  its  asymptotic  value, 
where  x_  — •  1/2.  More  precisely,  solving  the  quadratic  x  *  -  a(l  -  2x .).  where  a  is  given 

by  the  right-hand  side  of  (A- 22),  we  find  that  for  a  >  1  we  have  x_  ==  1  2.  In  particu- 

ji) 

lar,  for  a  =  3,  x_— 0  45.  Expressing  T  in  electron  volts,  letting  —  19  (llj  p  212), 

8 

Patm  ~  and  taking  the  logarithm  of  the  right-hand  side  of  (A-2‘2),  this  translates  to 


In ( 6.7  x  103)  -  ln(  1 .9)  -t  -2-ln(  1.161  T„v) 
Then,  letting  3  T,v  l„v.  we  find 


It 

T 


>  ln(3) 

PV 


(A-21) 


.^•-^-ln^  •  8.72+  -^-1  n I  (A-25) 

A  sufficient  condition  is  3-  5,  but  taking  the  ionization  potential  for  air  to  be 
1,  14  1  eV  (|i;  p  212),  we  find  3-  9.7.  Thus,  a  sufficient  condition  for  the  complete 
single  ionization  of  air  is 

T,v  >  (I,)ev/9-7  >  1.5  eV  (A-26) 

Note,  however,  if  the  temperature  rises  sufficiently  (  2  to  3  eV  ).  double  ionization  will 
commence  Thus,  within  the  temperature  range  1.5  eV  to  2  5  eV.  and  for  fre<|uencies 
satisfying  (tiuz),.,  ■  1.5  eV  (otherwise  wc  must  consider  bound-free  transitions),  we 

may  use  the  formulae  ( A- 19)  with  Z  --  1. 
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INCLUSION  OF  BOUND-FREE  TRANSITIONS 


Formulae  for  the  absorption  coefficient  which  include  the  effects  of  both  free-free 
and  bound-free  transitions  are  derived  semiclassically  in  |4i.  As  described  there  (  4  p 
266),  the  results  are  reasonably  accurate  in  a  more  general  setting.  We  note  that  in  their 


derivation  of  (5.34)  and  also  (5.41),  a  factor  of  2  — 


(rn+  1) 


t 


g 


|ml 


has  been  replaced  by  1  (a  rea¬ 


sonably  good  approximation)  We  reintroduce  this  factor  and  an  additional  factor, 
designated  £(<*.■),  which  is  simply  designed  to  indicate  a  quantum-mechanical  correction 
(When  free-free  transitions  dominate,  £(uu)  reduces  to  the  Gaunt  factor  g  See  1  .  p  241- 
242  )  Under  these  circumstances,  we  may  write  (|4j  p  278) 


G  -  E 


( A  -  2  7 


whet 


and 


•)  c . .  n...  _v 

(n>  df  — f  '  m  F  ,„(  x )  £,„(- 


T(io) 


T- 


( A-28a) 


16,t-  e" 
3v/3  lick2 


0  66  x  10  *  cur  deg" 


A- 28  b) 


m  Ini-1  kT;  x  li-'  kT 


(A-2SY 


x ) 


x  y 


F‘>r  x  v 


A-2S, 


Note  that  y.„  is  xtm  ol  4  The  ion  densities  are  determined  from  (A- 201  and  (A-21) 
The  absorption  coefficient  //,.  is  given  by  substituting  (A-27)  in  ( A  - 1 ) .  whnh  includes  the 
effects  of  stimulated  emission 

Several  observations  are  in  order:  First,  we  note  that  the  inequalities  in  (A-28<l) 
should  actually  read  and  “  >  >”  since  the  formulae  are  approximations  derived 

under  these  assumptions  [4j.  The  point  at  which  the  two  expressions  are  equal  is  actu¬ 
ally  x  ln(2  yni)  4  y,„,  but  this  occurs  quite  close  to  x  ym  Second,  an  assumption 
was  made  in  the  derivation  of  (A-28)  that  the  ions  may  be  modeled  as  hydrogen-like 
atoms  with  charges  of  Z;  this  may  not  always  be  accurate  (  4  p  275)  Also,  we  remark 
that,  in  the  derivation  of  (A-28),  [4]  sums  the  contributions  from  the  free-free  transitions 

t  This  really  should  be  the  ratio  of  the  ion  partition  unctions  u<ml,  blit  it  h.is  been  approximated 
bv  the  i.itio  of  the  statistical  weights  of  the  ground  states 
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and  the  bound-free  transitions.  Although  the  expression  for  the  free-free  transitions  is 
only  valid  for  tw  «  kT,  its  contribution  for  1w  >  kT  is  negligible.  Thus,  (A-28) 
remains  a  valid  approximation  under  quite  general  conditions;  i.e.,  so  long  as  <*/  »  wp 
and  w  »  f'm 

If  we  consider  the  case  of  single-level  ionization  with  m  =  0  (i.e.,  kT  «  Ij)  and 
frequencies  for  which  hcc  «  I,,  then  (A-28)  reduces  to  the  Kramers-Unsold  formula, 


16tr2  e6k  e  ,  ,  0  g'  -u  -  - 

3x/3  he  e°H"g(0)  (tiw)3e 


g|0)  (M 


Observe  that  the  conditions  for  the  validity  of  (A-29)  may  be  written 


(1  _e-1>WkT) 

(A-29a) 

r  (1  -  e“w/kT)  . 

(A-29b) 

mav  be  written 

(a)  kT  <  <  I, 

(b)  h  <jJ  <  <  I ] , 

which  is  less  stringent  that  (i').  Since  Z  =-  tn  +  1  =  1  and  n,  n_  —  n  2,  we  have 

_  P-2  no 

n°  _  (kT)2  nin_ 

2  2 

X.  p  n0 

(kT)2  n j n 


(  sing  this  expression  and  substituting  (A- 20)  into  (A-29),  we  find 


IL 


lf)7r*’es 


3\/3  hek  I  km 


3/2 


p-  1 


,.t>-  kT 


T‘ J* 


OH  f„(*1 


x  ”  Patm 


(T  I04)5  2  (I. 


~  (<■“• >T  i: 


1  ^  If- V 


(A-:K)a) 


(A- 30b) 


with  x,  computed  via  (A-22). 

Note  that,  up  to  the  edge  correction  e*1  (A-30b)  matches  equation  (6.46)  of  [1;. 

Also,  for  tiv*’  kT  and  x  1/2.  (A-30b)  reduces  to  (A-19d)  with  Z  1.  Finally,  we 
remark  that  in  comparison  with  (A-29b),  (6.15)  of  ;i!  seems  to  be  lacking  a  factor  of  two. 
This  occurs  even  if  we  set  the  factor  £0(^’)  ("hicli  for  our  purposes  is  actually  close  to 
two)  equal  to  one 
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